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We present a two-channel model to describe the quantum state of two atoms with finite-range 
interaction near a Feshbach resonance. This model provides a simple picture to analytically derive 
the wave function and the binding energy of the molecular bound state. The results agree excellently 
with the measurements and multichannel calculations. For small binding energies, the system enters 
a threshold regime in which the Feshbach molecules are identical to long range atom pairs in single 
channel. According to their threshold behavior, we find Feshbach resonances can be classified into 
two types. 
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I. INTRODUCTION 

Formation and Bose-Einstein condensation (BEC) of 
molecules Q have recently been achieved based on ultra- 
cold atoms with magnetically-tuned Feshbach resonances 
0- In these experiments, Feshbach coupling is induced 
by tuning a foreign molecular state near the scattering 
continuum, which allows for an efficient transfer of collid- 
ing atoms into molecules. This method works for virtu- 
ally all alkali atoms, and can create ultracold molecules 
from various sources including Bosc condensates yfl, de- 
generate Fermi gases Q , or normal thermal gases [3, El • 

Feshbach molecules have special and unique properties. 
They typically populate only one weakly-bound quantum 
state, and the bound state can strongly couple to the 
scattering continuum via Feshbach resonance. We may 
ask the following question: should Feshbach molecules 
rather be considered as molecules in a specific rovibra- 
tional state or as pairs of scattering atoms near the con- 
tinuum? This distinction is particularly crucial in the 
studies of the BEC to BCS (Bardeen-Cooper-Schricffcr 
state) crossover in degenerate Fermi gases, which call for 
a clarification of the quantum nature of the Feshbach 
molecules Q. 

Molecular states near Feshbach resonances have been 
recently investigated based on sophisticated and com- 
plete two-body or many-body theory 0, Q and multi- 
channel scattering calculations 0, El- All works sug- 
gest that the Feshbach molecule is generally a coherent 
mixture of the foreign molecule in the closed channel 
and long-range atom pair in the open scattering chan- 
nel. Near resonances with large resonance widths, the 
molecules can be well approximated as pairs in the open 
channel. For narrow resonances, as suggested by numer- 
ical calculation 0, E3 > the closed channel dominates and 
a short-range molecule picture is appropriate. 

In this paper, we use a simple two-channel model to de- 
scribe two interacting atoms near a Feshbach resonance 
(Sec. II). To account for the finite interaction range of 
real atoms, we introduce a spherical box potential, which 
allows us to analytically calculate the molecular bound 
state in different regimes and their threshold behavior 



(Sec. Ill and Sec. IV). Finally, we apply our model to Fes- 
hbach molecules in recent Fermi gas experiments and to 
characterize the associated Feshbach resonances (Sec. V). 



II. MODEL 

We model the interaction of two identical, ultracold 
atoms with mass m based on an open channel |o) that 
supports the scattering continuum and a closed chan- 
nel |c) that supports the foreign bound state. The 
wave function of the atoms is generally expressed as 
= V'oMlo) + Tp c (r)\c), where ip (f) and ipc{r) are 
the amplitudes in the open and closed channels, respec- 
tively, and r is the inter-atomic separation. We assume 
the interaction (h 2 /m)v is described by a spherical box 
potential with an interaction range of ro, see Fig. (1). For 
r > 7'o, the potential energy of the open channel is and 
the closed channel oo. For r < ro, the open (closed) chan- 
nel has an attractive potential of —h 2 q 2 /m (—h 2 q 2 /m), 
and a coupling term = |0| between the channels. The 
wave function satisfies the Schrodinger equation: 
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The solution of the above equation for zero scattering 
energy E = can be expressed as: 



for r > ro '■ 
for r < ro '■ \ip) oc 



r-a 

oc o 

r 

sin q + r 



A sin g_ r 



(3) 

R, (4) 



where the scattering length a and A are constants, q± 
are the "eigen wave numbers" for r < ro associated with 
the eigen states 
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FIG. 1: A two-channel model for Feshbach resonances. We 
assume a spherical box potential between two atoms with an 
interaction range of ro. A bound state with energy E c rela- 
tive to the scattering continuum is supported by the closed 
channel. 



III. SCATTERING LENGTH AND THE 
MOLECULAR EIGEN STATE 

Experimentally, the relative energy between the con- 
tinuum and the bare state can be adjusted linearly by 
a magnetic field B-induced Zeeman shift —fiB, where 
fj, = /jl — fi c and ii (ii c ) is the magnetic moment of the 
open(closed) channel. Replacing E c by E c + /iB, we can 
rewrite Eq. (8) in terms of the magnetic field as 
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where the resonance width AB and the resonance posi- 
tion B res are given by 



|+) = cos6>|o) + sin6>|c) 
|-) = -sin0|o) + cos6»|c), 

and tan 26 = 2VL/(q 2 - q 2 c ). 

Based on the boundary conditions ipc(fo) 
fpo(ro)/fp' { r o) = r Q (r /a - 1), we get 
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The latter equation shows how in general, each channel 
contributes to the scattering length. 

In cold atom systems, Feshbach resonances are, in most 
cases, induced by hyperfine interactions or spin-spin in- 
teractions. Both interactions are many orders of mag- 
nitude weaker than the relevant short range exchange 
potential. It is an excellent approximation to assume 
& < 9o> 1c an d \Qo — Qc\- Hence, we have 9 -C 1, q + ks q 
and g_ ~ q c . 

In this limit, the closed channel contribution is signif- 
icant only when the foreign state is close to the con- 
tinuum, in which case the last term in Eq. (7) di- 
verges. Given the energy of the closed channel state as 
E c = (h 2 /m)e c and e c <C Q c /ro, the boundary condition 
Vv (ro) = sin \J q 2 + f-c ro = allows us to expand the last 
term in Eq. (7) as — 7/e c . Here 7 = 2q 2 8 2 /r character- 
izes the Feshbach coupling strength. To the same order 
of expansion, the middle term in Eq. (7) is a constant 
across the resonance and can be identified as {tq — a^g) , 
where a bg is the background scattering length. Equation 
(7) reduces to 
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Several interesting features are shown here. First of 
all, we find the resonance width is proportional to both 
the Feshbach coupling 7 and the background scattering 
properties a\, g — ro- The latter dependence is due to the 
fact that the scattering amplitude at short range is pro- 
portional to the scattering length. A larger short range 
scattering amplitude leads to a stronger coupling to the 
closed channel. 

Secondly and importantly, the resonance position is 
offset by exactly AB relative to the crossing of the bare 
state and the continuum, B = —/.i~ 1 E c , sec Eq. (11). For 
a positive scattering length at g > ro , this shift is negative 
AB < 0. This feature leads to the "rcnormalization" of 
the Feshbach resonance location discussed in Ref. Q ■ 

To understand the origin of the resonance shift, we 
should return to Eq. (8). The divergence of the scat- 
tering length occurs when the open channel contribution 
(middle term) is exactly canceled by the closed channel 
one (last term). For systems with large background scat- 
tering lengths I a;, g I 3> ?*o and strong Feshbach couplings 
7, this cancelation can occur even when the bare state 
is far away from the continuum. A large resonance shift 
then results. 

Now we turn to the binding energy of the molecules. 
Assuming a bound eigen state \ip m ) exists near the con- 
tinuum at E = —E rn = —(ft 2 /rn)e m , where E m > 
is the binding energy, we can determine E m by follow- 
ing essentially the same calculation as Eq. (l)-(7). The 
equivalence of Eq. (7) gives 
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Assuming 8 <C 1 and the bound states in both channels 
are close to the continuum, namely, \abg\ 3> r o an d |e c | *C 
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nel, we can write the wave function of the eigen state as 



bare state detuning £ c (I ' ) 

FIG. 2: Generic energy structure of Feshbach molecules. 
Based on Eq. (13), we show the reduced energies of the con- 
tinuum (dotted lines), bare molecular states (dashed lines), 
and the dressed molecular state (e m , solid lines). We assume 
(a) cibg — ro — I and 7 = l~ 3 and (b) a,b g — ro = —SI and 
7 = i -3 , where I 3> ro can be any relevant length scale. Ar- 
rows mark the offset resonance positions. Insets show the 
threshold behavior of the bound state. 



3o/ r o> we can expand the two terms on the right side of 
Eq. (12) to leading order as (r — cifcg) -1 and — j(e c + 
£m)~\ respectively. Equation (12) then reduces to 



(e m + e c )( 
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This result shows the evolution of the eigen state near 
the resonance. Similar result is obtained in Ref. Q based 
on a contact potential. We can immediately see that in 
the absence of the Feshbach coupling 7 = 0, the solutions 
of Eq. (13) are e m = — e c and e TO = (at, g — ro)~ 2 (for 
a-bg > ro), which exactly correspond to the bare bound 
states in the closed channel and the open channel (for 
dbg > ro), respectively. 

In the presence of the Feshbach coupling 7 > 0, 
Eq. (13) suggests an "avoided level crossing" -like energy 
structure, see Fig. (2), which also illustrates the reso- 
nance position shifts. The level crossing, however, is not 
hyperbolic as it is in a two-level system. In particular, 
at small binding energies, the bound state energy ap- 
proaches the continuum quadratically, see Fig. (2) inset. 
Far below the continuum, the bound state approaches 
the bare state in the closed channel. 

To better quantify the role of the open and closed chan- 
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and the mixing angle (f> is defined below. 

We show in Eq. (15) that the eigen state generally 
occupies both the closed channel and open channel. We 
can introduce a mixing amplitude sin <p as the amplitude 
in the closed channel 
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The mixing fraction sin <f> can be evaluated by a di- 
rect integration of the closed channel wave function. Al- 
ternatively, noticing that the mixing also leads to a de- 
pendence of the eigen state on the bare state, we can 
also derive sin 2 from the dependence of E m on E c , 
or from the magnetic moment of the Feshbach molecule 
fi m = dE m /dB. All methods lead to the same result 
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THRESHOLD REGIME IN OPEN 
CHANNEL 
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Despite the seemingly complex equations shown in 
previous sections, the Feshbach molecules are simple 
and universal near the scattering continuum. Expand- 
ing Eq. (13) with small E m and using Eq. (9)-(ll), we 
find the binding energy of the Feshbach molecules has 
a simple dependence on the scattering length and in- 
creases quadratically in magnetic field near the reso- 
nance, namely, 
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where E* = {h/Am)^ 2 (a bg - r ) 4 . 

Equation (20) shows identical dependence on scatter- 
ing length and interaction range as of single channel 
molecules in the threshold regime [l2^ . Furthermore, 
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TABLE I: Parameters of the 6 Li and 40 K Feshbach reso- 
nances. Interaction range ro is derived from Ref. [T3 [ . see 
text. Feshbach coupling 7 is derived from Eq. (10). ao is 
Bohr radius and [ib is Bohr magneton. 
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taking the limit e m — > in Eq. (14) and (19), we find 
the molecular wave function here is purely in the open 
channel. Its spatial extent is determined by quantum 
uncertainty, (r) \J mE m — h/2, and can be much larger 
than the interaction range (r) = a/2 3> ro. In this limit, 
the Feshbach molecules are identical to long-range atom 
pairs in a single open channel. 

By expanding Eq. (19) at small e m and using Eq. (13), 
we find the closed channel fraction can be expressed as 



- - (1 + 8) — 

E* K ' E* 



(22) 



where S = j{a,b g — tq) 3 . 

From Eq. (22), we see that E* provides the leading 
order estimation of the closed channel admixture. When 
E m <C E* or E m <C E*/\S\ (this condition applies when 
S <C —1), the Feshbach molecule is purely in the open 
channel. As expected, the threshold regime is wider for 
resonances with larger 7 and \ab g — ro\- 

We can further determine the "open channel- 
dominated" regime by setting sin 2 </> < ^ in Eq. (19). 
For resonances with small |<5| <C 1, this condition corre- 
sponds to E m < E*, which, in terms of magnetic field, 
maps to only a small fraction of |<5| AB near the resonance 
B res . For resonances with large \S\ 3> 1, the open chan- 
nel dominates when E m < (7/2) 2 / 3 , which covers the full 
resonance width when ab g — ro < 0, and covers the entire 
upper branch of the bound state when a& s — 7"n > 0. 

Based on the range of the single channel regime, we 
suggest the broad(narrow) resonances be defined as those 
with \S\ > 1 (\S\ < 1). Within the width of the Fesh- 
bach resonance, the molecules associated with a broad 
(narrow) resonance are better described as long range 
pairs in the open channel (short range molecules in the 
closed channel). We note that this definition is purely 
based on two-body physics. 



V. FESHBACH MOLECULES IN b LI AND 4U K 

Finally, we apply our model to the 6 Li2 and 40 Fv2 Fesh- 
bach molecules created in recent BEC-BCS experiments 
0. These molecules are stable near the resonance and 
both the molecular binding energies and the scattering 
lengths have been well measured and studied [Til Il4j . 

To model the interaction of atoms, we adopt B res , AB, 
Bb g and \x from recent measurements and numerical cal- 
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FIG. 3: Binding energies £" m (dotted lines) and mixing frac- 
tions sin 2 (j> (solid lines) of the molecules near the 6 Li and 
40 K Feshbach resonances. The curves are calculated from 
Eq. (13), Eq. (19) and the parameters in Table 1. Binding en- 
ergies from multi-channel calculation [iH (dashed lines), from 
JILA group measurement \Tty (open square) and the mixing 
fractions measurement from Rice group [T6j[ (open circles) are 
shown for comparison. The shaded areas indicate the typical 
BEC-BCS crossover regimes, \a\ > 3000ao. 



culations. To account for the finite range of the atomic 
interaction, which at low temperatures is determined by 
the van der Waals potential of ~ — r -6 , we choose the 
interaction range ro in our model to be the mean scat- 
tering length a defined in Ref. 01 • This choice ensures 
the same behavior of the scattering length in the thresh- 
old regime |l2j| . All parameters arc given in Table I. 

In Fig. (3), we show the calculated binding energy E m 
and the mixing fraction sin 2 <fi of the Feshbach molecules 
for the two 6 Li resonances and one 40 K resonance. The 
results agree very well with the multi-channel calculation 
H3| and the measurements on molecular binding e nerg y 



:ncrg> 

m 



|lll [14J , magnetic moment |fj and mixing fraction 

Both the Li resonance at 834G and the K resonance are 
broad with S = —2900 and 4.7, respectively. The open 
channel-dominated regimes of E m < E*/\6\ = h 210 MHz 
for the 834 G Li resonance is also larger than the Fermi 
energy of Ep f=a h 20 kHz in the experiments. (Here the 
Fermi wave number is kp ~ (3000ao) _1 .) For the K res- 
onance, the full upper branch of the molecular state is 
open channel dominated with mixing fractions less than 
10%. Therefore, we conclude the open channel descrip- 
tion of these Feshbach molecules in the crossover regime 
to be a good approximation. 
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For the narrower Li resonance at ^543 G, we obtain 
S ~ 0.0005 and E* = h 31 Hz < E F = h 20 kHz. This 
indicates an extremely narrow open channel regime of less 
than 50 [iG near the resonance, where the gas parameter 
is still over kpa = 20. Crossover experiments based on 
these Feshbach molecules cannot be described by open 
channel atom pairs and may lead to qualitatively different 
physics. We attribute the large difference between the 
two Li resonances to their different couplings 7 and very 
different background scattering length a^g, see Table I. 

In the above discussions, we note that Fermi energy Ep 
is an external parameter which depends on the density 
of the sample. Whether the molecules in the crossover 
regime can be described by single channel strongly de- 
pends on the density. The 5 parameter, however, pro- 
vides a better and independent measure to classify Fes- 
hbach resonances. We find that the two Feshbach reso- 
nances in 6 Li are the two extremes of broad and narrow 
resonances with S = —2900 and 5 = 0.0005. 



In summary, the two-channel model provides a simple 
picture to understand the molecular state near the Fes- 
hbach resonances. The analytic results of the molecular 
binding energy and mixing fraction on 6 Li and 40 K agree 
with the measurements and other sophisticated calcula- 
tions very well. Based on the threshold behavior of the 
bound state, we suggest a dimcnsionless parameter to 
assess the "broadness" of the Feshbach resonance. 
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